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We study (3+l)-dimensional holographic superconductors in Einstein-Gauss-Bonnet gravity both 
numerically and analytically. It is found that higher curvature corrections make condensation harder. 
We give an analytic proof of this result, and directly demonstrate an analytic approximation method 
that explains the qualitative features of superconductors as well as giving quantitatively good nu- 
merical results. We also calculate conductivity and u> g /T c , for uj g and T c the gap in the frequency 
dependent conductivity and the critical temperature respectively. It turns out that the 'univer- 
sal' behaviour of conductivity, ui g /T c ~ 8, is not stable to the higher curvature corrections. In 

' the appendix, for completeness, we show our analytic method can also explain (2+l)-dimensional 

' superconductors. 
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C/5 ■ I. INTRODUCTION 

It is often felt that the most remarkable discovery in string theory has been the AdS/CFT correspondence, [l[, which 
has been further extended to the gauge/gravity correspondence [2|- Interestingly, the gauge/gravity correspondence 
may play an important role in condensed matter physics 0, 0|- I n particular, the application of the gauge /gravity 
correspondence to superconductors has been intensively studied 0, @, 0, H, 0, E3, EL [HI ED, 0, [H, OS EH 
(see recent lecture notes [1, El for complete references). It would be very exciting indeed if we could explain high 
' temperature superconductivity from black hole physics. In addition, from the gravity perspective the existence of 
, scalar condensation in black hole systems itself deserves further stud y in relation to the "no-hair" theorems and a 
better understanding of the dressing of horizons by quantum fields p, Il9l|. 

Remarkably, on the gauge theory side, there is a puzzle. As is well known, the Mermin- Wagner theorem forbids 
continuous symmetry breaking in (2+l)-dimensions because of large fluctuations in lower dimensions. Nevertheless, 
£C) I holographic superconductors are found in (2+l)-dimensions. It is possible that fluctuations in holographic supercon- 
. ductors are suppressed because classical gravity corresponds to the large N limit. If this is true, then higher curvature 
corrections should suppress condensation. Of course, to examine whether or not the Mermin- Wagner theorem holds, 
we need to study 4-dimensional higher curvature gravity. Unfortunately, higher curvature gravity in 4 dimensions 
t^- ' is not particularly illuminating: higher derivative terms in general introduce ghost degrees of freedom [2(| . the ex- 
ceptions being either Gauss-Bonnet or Lovelock gravity [2l|, in which specific combinations of the curvature tensors 
are used, or f(R) gravity, [22], in which powers of the Ricci scalar only are used. Unfortunately, the former case is 
non-dynamical in 4 dimensions, and the latter case is conformally equivalent to scalar-tensor gravity, (23[, and black 
hole solutions are therefore identical to the Einstein case [24[ ■ 

To explore this issue, we instead study 5-dimcnsional Einstein-Gauss-Bonnct gravity, which gives a known gen- 
eralization to the Schwarzschild black hole solution (2f|. We would also like to investigate if the universal relation 
between the gap uj g in the frequency dependent conductivity and the critical temperature T c : uo g /T c ~ 8, found in 

0, is stable under stringy corrections. In the case of the quark-gluon plasma, there is a universal shear viscosity 
to entro py density ratio rj/s = \/Air [27]], and there are several analyses investigating the stability of this universal 
relation [28l. l29l. l3fj|. l3ll. [32. l33l. 13-iI [3fjj to higher curvature corrections. To the best of our knowledge, no corresponding 
analysis exists in the case of superconductors. Hence, we look at gap frequency at a given temperature numerically 
to explore its stability under higher curvature corrections. 

To investigate the effect of the higher curvature corrections on the superconductor, we operate in the 'probe' limit, 

1. e. where the gravitational back reaction of the scalar and vector fields on the background geometry is neglected. 
At least for temperatures near the phase transition this should be a good approximation, and has been found to be 
valuable in the Einstein limit Q . Ideally, one would like to have a full analytic description of the phase transition and 
condensation phenomena, and in this paper we take a modest first step in this direction. We first prove the existence of 
a bound on black hole temperature above which no condensation can occur. Since there is always an analytic solution 
with vanishing scalar, Q, we cannot similarly prove the existence of a nontrivial scalar solution below T c , however, a 
simple matching method provides an approximate analytic solution which explains the phase transition behaviour and 
gives a very good approximation to the phase diagram. Indeed, we can calculate the critical temperature analytically 
within a few percent in the best case. In a sense, this is the most important result in our paper. Numerical methods 
complete the proof of condensation, and are clearly necessary for fully describing the properties of the fields and the 
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details of the physics. 

The organization of the paper is as follows. In section II, we introduce the model and numerically demonstrate the 
effect of the Gauss-Bonnet term on the superconductor. We find that stringy corrections make condensation harder. 
In section III, we present an analytic explanation of the superconductor. We can understand the qualitative features 
of the superconductor with a simple calculation. The analysis also gives fairly good numerical results. In section IV, 
we study the conductivity and show the universality is unstable under the stringy corrections. We conclude in section 
V. In the appendix, we present an analytic explanation of (2+l)-dimcnsional superconductors for completeness. 



II. GAUSS-BONNET SUPERCONDUCTORS 



In this section, we study the effect of Gauss-Bonnet term on the (3+l)-dimensional superconductor using the probe 
limit. In the probe limit, gravity and matter decouple and the system reduces to the Maxwell field and the charged 
scalar field in the neutral black hole background. 

We begin with the Einstein-Gauss-Bonnet action: 
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(1) 



where g is the determinant of a metric g^v and R^vXp, R^v and R are the Riemann curvature tensor, Ricci tensor, 
and the Ricci scalar, respectively. We take the Gauss-Bonnet coupling constant a to be positive. Here, the negative 
cosmological constant term —6/L 2 is also introduced. The background solution we consider is a neutral black hole [25| : 



ds 2 



-f{r)dt 2 



—— + —r{dx 2 + dy 2 + dz 2 ) 
f(r) L 1 
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where 
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Here, M is a constant of integration related to the "ADM" mass of the black hole j26|. The position of the horizon 
defined by f(rn) = is at rjj = {ML 2 ) 1 / 4 ". In order to avoid a naked singularity, we need to restrict the parameter 
range as a < L 2 /4. Note that in the Einstein limit (a — ► 0), the solution goes to f(r) = — 72-, and L can be 
regarded as the curvature radius of asymptotic AdS region (r — > 00). For general a, however, the solution Q behaves 
as 
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in the asymptotic region. Hence, we define the effective asymptotic AdS scale by 
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The Hawking temperature is given by 
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where a prime denotes derivative with respect to r. This will be interpreted as the temperature of the CFT. 
In this background, we now consider a Maxwell field and a charged complex scalar field, with the action 



S = / d 5 xv^g~ 



--F^F^ - |W - iAi\>\ 2 - m 2 |^| 2 
Taking a static ansatz, = (4>{r), 0, 0, 0, 0) and ijj = ijj{ r )> the equation of motion for (f>(r) becomes 

3 2y, 2 

(p + -<P = 
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where without loss of generality ip can be taken to be real, and satisfies 



Note that the Maxwell equations imply that the phase of i/> must be constant, which is set to zero by a residual gauge 
for A^. 

We now want to solve (J3J) and © for the scalar and vector field. For the main part of this paper, we choose to set 
the mass of the scalar field to be m 2 = — 3/L 2 , so that the mass remains the same as we vary a. Note however, that 
because of the variation of the effective asymptotic AdS curvature, with a relative to L means that this mass 
actually increases (i.e. becomes less negative) with respect to the asymptotic AdS scale. On the other hand, while 
setting to 2 = — 3/i 2 ff has the advantage of fixing the mass relative to the asymptotic AdS scale, this mass now varies 
with respect to the physical measurablcs of black hole mass and temperature as we vary a. Since condensation is 
a temperature dependent phenomenon, we believe that fixing the scalar mass with respect to the black hole is the 
correct physical choice, however, we have also checked that for the alternative choice of mass the same qualitative 
features occur as we vary a. 

In order to solve our equations we need to impose regularity at the horizon and the AdS boundary: 
• Regularity at the horizon gives two conditions: 

<Kr B ) = 0, il){r H ) = ~r H ^'(r H ) . (10) 



Asymptotically (r — *■ oo) the solutions are found to be: 



where A± = 2± \J A — 3 (■^ff 1 ) 2 - Here, p and p are interpreted as a chemical potential and charge density, respectively. 
Note that these are not entirely free parameters, as there is a scaling degree of freedom in the equations of motion. 
As in [||, we impose that p is fixed, which determines the scale of this system. For if), both of these falloffs are 
normalizable, so we can impose the condition either C_ or C+ vanish. We take C_ = 0, for simplicity. 

According to the AdS/CFT correspondence, we can interpret (O) = C+, where O is the operator dual to the scalar 
field. Thus, we are going to calculate the condensate (O) for fixed charge density. The results are shown in Figure 
[TJ From Fig. [TJ we see the GB term makes the condensation gap larger. We also see that the Chern-Simons limit 
shows a slightly different dependence of the condensate on temperature. This can be understood from the behaviour 
of gravity near the horizon. In the Chern-Simons limit, a = L 2 /4, we get 



„, . 2r 2 ( VML 2 \ , . 

= T? y r 2 J (12) 

Hence, the correction to the AdS quadratic gravitational potential dependence is simply a constant instead of a 1/r 2 
dependence, leading to more gentle tidal behaviour. The process of scalar condensation (or the formation of scalar 
hair) can be understood as arising in part from the 'negative' mass of the scalar field, but also as arising from the 
potential well that occurs near the horizon. For black holes with large mass, this well is too broad and shallow to 
allow for the formation of a nonzero scalar, however, for small black hole mass, the strong curvature near the horizon 
is amenable to condensation. (Refer to the analytic arguments in section III which show how the behaviour of the 
gravitational potential interacts with features of the scalar condensate.) At some stage further decreasing the mass 
of the black hole does not alter the shape of the condensate much, as the scalar is already sampling regions of strong 
curvature. However, the CS limit has a rather different and smoother profile near the horizon, therefore it is not 
surprising that decreasing the black hole mass in this case has more impact on the details of the scalar field. 

Numerically, we found that increasing a resulted in a decrease of the critical temperature: T c = 0.198/? 1 / 3 for 
a = 0.0001, T c = 0.186/9 1 / 3 for a = 0.1, T c = O.tflp 1 / 3 for a = 0.2, and T c = 0.158/9 1 / 3 for a = 0.25 (see also figure 
[5]). Thus the effect of a is to make it harder for scalar hair to form. Changing the scalar mass to m 2 = — 3/L 2 ff 
gives a similar, though less marked, behaviour, for example T c = 0.181/9 1 / 3 for a = 0.2. We can therefore conclude, 
as expected, that the higher curvature corrections make it harder for the scalar hair to form. One can expect this 
tendency to be the same even in (2+l)-dimensions, however, it remains obscure to what extent this suppression affects 
the physics of holographic superconductors in (2+l)-dimensions. 

We have thus numerically verified that Gauss-Bonnet superconductors exist. However, we would ideally like to have 
an analytic understanding of condensation to back up this numerical work. This is what we now turn to. 
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FIG. 1: The condensate as a function of temperature for various values of a. The (lowest) red line is for a — 0.0001, the 
middle brown plot is ot = 0.1, the top blue line is a = 0.2 and the remaining line intersecting the other three in green is the 
Chern-Simons limit a — 0.25. Note that while the generic Einstein-Gauss-Bonnet behaviour is to level out for T < T c /2, the 
Chern-Simons limit has a much stronger variation of the condensate with temperature. 



III. SUPERCONDUCTORS IN A NUTSHELL 



Although in the previous section we used numerical integration to explicitly demonstrate the condensation phe- 
nomenon, ideally we would like to obtain an analytic understanding in parallel. Since our equations are nonlinear 
and coupled, we cannot derive analytic solutions in closed form, however we can deduce a great deal of information 
analytically. We first prove the nonexistence of condensation for large T before explicitly deriving the phase diagram 
analytically by using approximate solutions. 

Note that the trivial solution to (jSJ and 



WD = 4(l-J) 



(13) 



ip = (14) 

always exists. We will now prove that there is a temperature above which this is the only solution. 

First consider the <f) equation (J5J). Let <p(r) = <fro(r) + 8(f), where (f>o(r) is defined above. Then (JSJ) implies 

(r 3 ^')' > (15) 

however, as r — ► oo, r 3 <p' — > 2p = r 3 <fi' , hence r 3 8<fi' — > at infinity, and using 5(f) = at rji we have that 8(f) 1 < 0. 
Hence 

0(r) < Mr) ■ (16) 
Next consider the scalar field, and define the variable X = rip: 

Now, the boundary conditions at the horizon imply X' H = Xn/ArH, and at infinity, rfX' — > 0, thus the existence of 
a condensate requires a turning point in X, X'(r T ) = 0, with X" < for X > 0. This in turn requires 

4>l{r T ) 3 f(r T ) f(r T ) 2 (r r ) 3 f(r T ) f(r T ) 

f(r T ) L 2 r T r\ > f(r T ) + r T r% 1 ° J 

at the turning point. By inputting the form of (f>o(r), it is easy to see that if M is too large, this inequality can never 
be satisfied, as the combination of (f> and the geometry to the LHS of (fTS)) is always negative. This gives a loose upper 
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FIG. 2: A comparison of analytic and numerical results. The shaded region is that in which the geometry forbids the possibility 
of a scalar condensate from (|18|) . The dashed line indicates the analytic approximation of the value of T c obtained by matching 
methods, (|39p . The data points are the exact numerical results. For simplicity p and L have been set to 1. 



bound on the critical temperature as shown in figure [2] (For a = we need to use the fact that J (rfX') = to 

bound M. This also gives a tighter analytical bound for nonzero a, however, the above argument is more direct.) 

We have thus numerically verified that Gauss-Bonnet superconductors exist. Having shown that there is a critical 
temperature below which there is no barrier to condensation, we will now show we can understand the essential 
features of condensation by using approximation techniques. 

Once again, let us change variables and set z = Under this transformation equations ([8]) and ([9]) become 

<t>"--/- J^0 = O (19) 

." + (H>+*(^w>=° 

where a prime now denotes 4-. The region rn < r < oo now corresponds to < z < 1. The boundary conditions 
now become: 

• Regularity at the horizon z = 1 gives 

0(1) = 0, = • (21) 

• In the asymptotic AdS region: z — -> 0, the solutions are 

<j> = /.i-qz 2 , ip = D„z x - +D+z x + , (22) 

where X± is the same as in equation pip . As boundary conditions, wc fix qr 2 H and take Z?_ to be zero. 

We now find leading order solutions near the horizon and asymptotically, say 1 > z > z m and z m > z > 0, and 
then match these smoothly at the intermediate point, z m . As a consequence, we will demonstrate the phase transition 
phenomenon directly, and derive an (approximate) analytic expression for the critical temperature. Moreover, we will 
have a much better analytical understanding of a dependence of the critical temperature, as the proof above only 
gives a loose bound on the critical temperature and only indirect access to an expression. 
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A. Solution near the horizon: z = 1 



We can expand and ip in a Taylor series near the horizon as: 

<P(z) = 0(l)-0'(l)(l-z) + ^"(l)(l-z) 2 + - 
1>{z) = V(l)-^(l)(l-^) + ^"(l)(l-^) 2 + 



(23) 
(24) 



From (|2Tj) . we have 0(1) = and ^>'(1) = jip(l), and without loss of generality we take </>'(!) < 0, "0(1) > to have 



(z) and 4>{z) positive. Expanding (fl^) near z = 1 gives 
= 0'(1) 



1 ,/ i ri 2^ 2 
2rf^(l) 2 



* 4 (1 -,)/'(!) I + + 



Thus, we get the approximate solution 

4>{z) = -0'(i)(i - z ) + \(i- y^(i) 2 ) 0'(i)(i - + 

Similarly, from (f20|) . the 2nd order coefficients of 4> can be calculated as 



V/'(l) = -V' 



0'(1) 



2 W + 39-^ 
^7 



zV 2 







2=1 



4z 3 /V + z 4 /'V + z 4 /V" + 3?$V' 



4z 3 / + z 4 f 



r^(-^(l)(l-z) + ...f 

/'(!)*(! 



z = l 



#'(1) 2 *(1) 



where we used THopital's rule at the second term in the second line. Thus, we get 
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(25) 



(26) 



(27) 



(28) 



8 I? 

After eliminating ^'(1) from above equation by using Eq. (|2ip . we find an approximate solution near the horizon as 

L 4 



^) = 4^ + 4^ + (-64 + 2^ 64,! 



*'(1) 2 U(1)(1 -*)* + • 



n 



(29) 



B. Solution near the asymptotic AdS region: z = 

From (|22p , and ip in the asymptotic region are given by 

<j>{z) = [i-qz 1 , 0(z) = D + z x+ 

where qr 2 H is fixed and we have set D_ = from the boundary condition. 



(30) 



C. Matching and Phase Transition 

Now we will match the solutions (|26|) . (|29|) and ((30|) at z m . Interestingly, allowing z r „ to be arbitrary does not change 
qualitative features of the analytic approximation, more importantly, it does not give a big difference in numerical 
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values, therefore for simplicity in demonstrating our argument we will take z n 
asymptotic solutions smoothly, we require the following 4 conditions: 
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1/2. In order to connect our two 

(31) 
(32) 

(33) 
(34) 



where we have set "0(1) = a an d — 4>'0-) = 6 (a, 6 > 0) for clarity. Now, the AdS/CFT dictionary gives a relation 
(O) 



LD + r^ L 2A +, hence we need to compute D + . From (|33j) and (p4|) we obtain 
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8 A+ + 2 



Using (|3"Tj) and (f52"| . a is expressed by 



where b is obtained from (f53"|) and 



Aq_ 

m 



b 
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assuming / (i.e. the scalar solution is non-trivial) as: 



6 = 8- 



3A, 
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(35) 



(36) 



(37) 



2(A+ + 2) 
e 1 

Using the Hawking temperature T = we can rewrite 

(38) 



Now we go back to the original variable, r, and compare the results with those in Q. First of all, we should note 
the relation p = q r%. We also define b by b = brn/L 2 . 
([55]) as 



L 2 T 3 



rp3 

'"if 



where we have defined T r as 



bL 



1/3 



1 



(39) 



We can now read off the expectation value (O) from (|35[) and ([55)) as: 



T, 



2tt 



13 
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(40) 



where we have normalized by the critical temperature to obtain a dimcnsionlcss quantity. We find that (O) is zero at 
T = T c , the critical point, and condensation occurs for T < T c . We also see a behaviour (O) oc (1 — T/Tc) 1 / 2 which 
is a typical mean field theory result for a the second order phase transition. 

Next, we evaluate the critical temperature from (|31)|) . The value of T c is 0.201p 1 / 3 /L when the Gauss-Bonnet term 
is absent, this should be compared with the numerical result T c = 0.198p 1 / 3 /L in Q- We therefore see that our 
analytic approximation is good. Moreover, as a increases to 0.1,0.2 and 0.25, T c decreases to 0.196,0.191 and 0.188 
respectively, which is in good agreement with our numerical results. 

Thus, we have (approximately) reproduced our numerical results from a simple analytic calculation. In particular, 
we have calculated extremely good estimates of the critical temperatures, and revealed how the structure of the 
interaction term has produced the phase transition. 
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IV. CONDUCTIVITY AND UNIVERSALITY 



We now calculate the conductivity, a, of our boundary theory. In [7|, the conductivity for various cases was 
calculated and it was found that there is a universal relation 

(41) 

with deviations of less than 8 %. The purpose of this section is to examine if this universality holds in the presence 
of stringy corrections. 

As A M in the bulk corresponds to the four-current J M on the CFT boundary, we can calculate the conductivity 
by considering perturbation of . The spatial components of A^ are decomposed into longitudinal and transverse 
modes: Ai = (diX,Aj-). These linearized perturbations are decoupled from each other and can be studied separately. 
The linearized equation of motion for Aj~(t, r, x 1 ) — j4(r)e ik x_iwt ei, which corresponds to the current density, is 

/' ,l\„. k 2 2,. a 



A " + K7 + r) A ' + \p~f-r) A=0 - m 

We solve this under the following boundary conditions near the horizon: 

A(r) ~ f(ry 1 ^ , (43) 

which corresponds to no outgoing radiation at the horizon. In the asymptotic AdS region (r — > oo), the general 
solution takes the form 



A = A (0) ^ AW^-tf^logAr 



(44) 



where A<°), A^ and A are arbitrary integration constants. Note the appearance of the arbitrary scale A, which leads 
to a logarithmic divergence in the Green's function, as explained in Q. Since this can be removed by an appropriate 
boundary counterterm, this scale will disappear from the results. 

From linear response theory, the conductivity can be calculated by the formula 

<t(w) = — G R (u>,k = 0) , (45) 
iuj 

where k is the wavenumber. The retarded Green function G R can be calculated through the AdS/CFT correspon- 
dence H(| as: 

G R = - lim f(r)rAA' . (46) 

r — >oo 

Thus, by using the solution (|44[) . the conductivity is given by 

k=0 Z 



We therefore need to solve (J42j numerically with the boundary condition (|43|) to obtain ^4^°^ and A^ asymptotically. 

The plots in figures 3—6 show the results of this numerical integration for a = 0.0001,0.1,0.2 and 0.25 at tem- 
peratures T/T c w 0.152,0.151,0.152 and 0.152, respectively. The red line represents the real part, and blue line the 
imaginary part of a. Taking look at the imaginary part of the conductivity, we see a pole exists at w = 0. From the 
Kramers-Kronig relations, this implies the real part of the conductivity contains a delta function. 

Clearly, the real part of the conductivity shows a frequency gap which indicates a gap in the spectrum of charged 
excitations. As a increases, the gap frequency (normalized by T c ) becomes large. As we noticed that condensation is 
an increasing function of a this tendency is consistent with the conventional relation tu g oc (Q) . We see the universal 
relation ^ ~ 8 found in Q is unstable in the presence of the Gauss-Bonnet correction. We have also checked that 
this conclusion is not affected by choosing the alternative scalar mass, M 2 = — 3/L 2 s . 
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FIG. 3: Conductivity for a = 0.0001 case 



FIG. 4: Conductivity for a — 0.1 case 





FIG. 5: Conductivity for a = 0.2 case 



FIG. 6: Conductivity for a = 0.25 case 



CONCLUSION 



We have studied holographic superconductors in the presence of Gauss-Bonnet corrections to the gravitational 
action. Motivated by the Mermin- Wagner theorem, we have investigated if the higher derivative corrections suppress 
the phase transition or not. We numerically solved the system in the probe limit and obtained phase diagrams for 
various Gauss-Bonnet couplings a, and calculated the critical temperatures. As we increase a, the critical temperature 
decreases, thus it turns out that stringy corrections make condensation harder. However, we did not reach the point 
that the critical temperature of the transition vanishes for changing a. We would also expect this to apply to the 
(2+l)-dimensional case, however, it is beyond the scope of this paper to determine if this could destroy holographic 
superconductors in (2+l)-dimensions. 

To understand phase transition phenomena, we also conducted an analytic analysis of the coupled nonlinear equa- 
tions, finding an approximate analytic solution. In spite of the apparent crudity of this approximation, we have 
analytically demonstrated the phase transition. Surprisingly, it turned out that the analytical method gave good 
agreement with the numerical results. In particular, we have calculated the critical temperature analytically. We 
obtained T c — 0.20T/9 1 / 3 , which is close to the numerical result T c = 0.198/3 1 / 3 for a — ► 0. We also applied the same 
method to the (2+l)-dimensional superconductor, presented in the appendix. The resultant critical temperature was 
T c = 0.103^, which should be compared with the numerical result T c = 0.118^/p 0]. 

Our other purpose was to examine the universality of the gap frequency to the critical temperature ratio. By 
calculating conductivity, we found the universal behaviour of conductivity lo 9 /T c ~ 8 was unstable to the stringy 
corrections. 

There are many issues to be investigated further. The obvious next step is to incorporate back reaction, which 
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is particularly important in the low temperature regime which corresponds to small black holes. In that case, the 
stability of black holes should be considered [37], Ha, Although we have investigated the stability of the 

superconductor under stringy corrections, it is also in trig uin g to study the dynamical stability of the condensation 
phase, as well as other aspects of superconductors fill |42| l43l |44| . 
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APPENDIX A: ANALYTICAL APPROACH TO (2+l)-DIMENSIONAL SUPERCONDUCTORS 

In the main text, we have shown a simple analytic treatment gives a good explanation of superconductivity. Here 
for completeness, we show (2+l)-dimensional superconductors can be also explained using the same method. 
In the 4-dimcnsional case, we have the following equations 

— 2 ^ = 0, (Al) 



/ 

/' 2\ (<? m 2 



^+(V + r)^+(^-^)^ = 0, (A2) 



where now 



with th = {ML 2 ) 1 / 3 . We set the mass of the scalar field, m 2 = — 2/L 2 , as in [f|. By changing to the z variable as 
before, z = r -f, (|A"Ijl and (|A"2]) become 

- ^-0 = (A4) 



z 2 {l-z 3 ) 

where a prime now denotes Next we consider the boundary conditions with these new variables. Regularity at 
the horizon, z = 1, requires 

0(1) = 0, ^'(l) = |v(l), (A6) 

and the asymptotic solution in the AdS region, z — > 0, reads 

<f> = fj,-qz, ip = C±z + C 2 z 2 . (A7) 

As in Q we fix the charge qrn and take C\ to be zero. 

We now find an approximate solution around both z = 1 and z = using Taylor expansion as before, then connect 
these solutions between z ~ 1 and z = 0. 



We expand </> and ip as 



Solution near the horizon: z — 1 



4>{z) = 0(1) - 0'(1)(1 -z) + ^"(1)(1 - zf + ■ ■ ■ (A8) 
*{z) = ^(1) - ~z) + ^"(1)(1 - z) 2 + • • • (A9) 
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From the boundary condition (|A6[) . 0(1) = and ip'(l) = ^tp(l), and we again set 4>'(1) < and V'(l) > for 
positivity of <f>(z) and ip{z). 

First we compute the 2 nd order coefficient <fi using (|A4[) as 



2L 2 ^ 2 
'* =1 ~ z 2 (l-z 3 )' 



■-iY(W>0> 



giving 



(z) = -0'(1)(1 -z)- -LV(1)V(1)(1 - ^ + 



The 2nd derivative of ^ is calculated similarly as 



z=l 



z(l-z 3 ) 

(z 4 + 2z)^" + 4z 3 V 
2z - 5z 4 



z 2 (l-z 3 ) 



L 4 



r=1 r\{l-z^ 



1> 



z=l 



L 4 



Thus 



^"(i) = -|^(i)-^r^( 1 )V(i)- 



Using (|A6|) to eliminate ip', we find 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



2. Solution near the asymptotic AdS region: z = 

We expand </> and i/j, making use of asymptotic solutions (|A7j) . as 



0(z) = n - qz + ^"(0)z 2 + 
1>{z) = C 2 z 2 + --- 



where we have used C\ = 0. 

Then the 2 nd derivative of cj) is given by 



2LV 



= 



2 = 



and we get simply 



(z)=fx — qz, ip(z) = C 2 z 2 



(A15) 
(A16) 



(A17) 



(A18) 



where qrn is fixed. 
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3. Matching and Phase Transition 

As before, we connect the solutions (|A11[) , (|A14[) and (|A18[) at z = i . In order to connect those solutions smoothly, 
we require the following 4 conditions: 

H- l -q= l -b+^a 2 b (A19) 

-q=-b- — a 2 b (A20) 



1 „ L 4 

"h 

C 2 = ^a + ^-ab 2 (A22) 



C 2 = —a ^-ab 2 (A21) 

4 18 lUr 2 



9 36r 2 H 



where "0(1) = a an d — 0'(1) = ^ with (a, 6 > 0) as before. Eliminating a 2 b from (|Al9j) and (|A20|) gives 



»=\q+\b. (A23) 



From (fA19l) and dA20l) . we can also deduce 



,2 



12 
Z25 



. (A24) 



The above relation alludes to phase transitions, namely, given q, fi has a maximum value when we assume the 
non-trivial solution a ^ 0. Substituting the relation (|A23[) into (|A24[) . we have 



To relate this result to the expectation value of the dimension 2 operator (O2) = \f2Cir 2 H jL z , we eliminate a6 2 from 
(TA"2~T]) and (fA"22j) to obtain 

C 2 = |o . (A26) 
Similarly, eliminating C2 from (|A21j) and ()A22|) gives 

- •>*=, 

L 4 



a I b 2 - 28-f 1=0, (A27) 



which determines b = 2v / 7rw/£ 2 provided a^O. 

Now we are in a position to reveal the phase transition phenomenon in this simple system. Noting the relation 
P = q th, and using the Hawking temperature: T = j^f?, (O2) can be expressed by 



where T r is defined as 



T c = 3 ^1_ (A29) 

We see that (O2) is zero at T = T c , which is a critical point, and condensation occurs at T < T c . The mean field theory 
result (0 2 ) oc (1 - T/T c ) 1 / 2 is also recovered. The value, (|A29[) . of T c is evaluated as 0.103y/p/L. Comparing with 
the numerical result 0.118-^p/L in Q, we find our analytic approximation is quantitatively good. Also the coefficient 
of (1 - T/T c y/ 2 as T -> T c is now 101T C 2 , while the numerical result is 144T 2 Q, which means this approximation 
seems good. 
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One may wonder what happens if we change z m . If we connect the solutions at z m (0 < z m < 1), the result is 



, , 16tt 2 2 + z m I 3 / T / Y , , 



where 



r.-JUf, S=.ffi±^ (A3i) 

4ttL \ b \j 1- z m 

In order to get the same value T c = 0. 118 y/p/L as @, we need to choose z m = 0.34. For this value, the coefficient of 
(1 — T/Tc) 1 / 2 as T — ► T c becomes 121T 2 . Thus, a numerically better approximation is possible, however, the choice 
of z m does not give a big qualitative difference. 
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